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Abstract 

. In the present paper problem of vortex images in annular domain between two coaxial cylinders is solved by 

' q-elementary functions. We show that all images are determined completely as poles of q-logarithmic function, 

where dimensionless parameter q = r\lr\ is determined by square ratio of the cylinder radiuses. The resulting 
solution for the complex potential is represented in terms of the Jackson q-exponential function. By composing 
pairs of q-exponents as the first Jacobi theta function and conformal mapping to rectangular domain we link 
our solution with result of Johnson and McDonald. We found that one vortex cannot remain at rest except at 
' ^ ' ■ the geometric mean distance, but must orbit the cylinders with constant angular velocity which is related with 

^ q-geometric series. Vortices in two particular geometries in g — > cxa limit are studied. 

^ ; 1 Introduction 

The classical method of images introduced by W. Thomson in 1845 becomes powerful method for solving boundary 
value problems in electrostatics and hydrodynamics [I],[2]- The method has been successfully applied to simple 
geometries as spheres, cylinders and half-spaces, where explicit formulas have been given. Unfortunately for 
complex body shapes the image principle becomes extremely difficult even to find approximate solution. This is 
why the image problem for which solution can be found in an exact form, like merging cylinders [3] for example, 
, becomes a member of very exclusive family . In the present paper we consider exact solution of planar vortex 
problem in annular domain between two coaxial cylinders by method of images. This problem has many interesting 
, applications. One of them is related to hydrodynamic interaction in which the modification of ambient fiow by 
' cylinders is carried out by obtaining the effect of single cylinder on the fiow and then applying the boundary 
, conditions on each cylinder to determine the unknown coefficients that appear in the series expansions. [1] solved 
' the diffraction problem of water waves by multiple cylinders placed at the free surface. 

Another application is related with inviscid two-dimensional fiuid dynamics experiments with magnetized elec- 
tron columns confined in a cylindrical trap[19j. The fiow vorticity is proportional to the electron density and the 
electric potential is analogous to the two dimensional stream function. Thus, the electrons mimic ideal two dimcn- 
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■ sional fiuid equations and by creating electron columns with the appropriate density, one can model fiuid fiows. It 
allows to model real problems of vortex interaction with topography [9] . Motion of a vortex in the neighborhood of 
a cosmic string [20.J and infiuency on this the cylindrically compactificated extra space dimensions is another class 
of possible cosmological applications. 

The problem of one vortex and one cylinder is connected with the Circle Theorem of Milne- Thomson [1] which 
can be rewritten for the complex velocity of the flow Viz) = ui ~ iu2 in the form 

= (1) 

where v{z) is complex velocity of the flow in unbounded domain, and the second term in ([T]) represents the correction 
to the complex velocity by the cylinder of radius ri placed at the origin. For a vortex at zq, of strength k and 
circulation, T = — 27rK, ([1]) can be written explicitly as 

V{z) ^— ^ + - (2) 

Z ~ Zq z — — Z 

zo 

where the second term represents a vortex of strength — k at the inverse point of zq, with respect to the cylinder. 

Zq 

Henceforth, we shall call the vortices at inverse points and at the centres of cylinders (or at the infinity) " vortex 
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images" or simply "images" . Therefore, in ([5]), there are two images; one positive image at the centre of the cyhnder 
and another negative image at the inverse point. In fact, images are used to replace the circle in the infinite 2-D 
plane. 

Another application is the case of a vortex at point zq inside a cylindrical domain with radius r2, C : \z\ < r2, 

V{^) = 2 (3) 

za 

where the vortex image is located at point rl/zo outside C. The solution ([3]) can be obtained from the circle theorem 
by first using the mapping z = 1/lo and also from the Laurent series expansion of the solution. 

The above two examples are limiting cases of the problem of a point vortex in annular domain between two 
coaxial cylinders with inner radius ri and outer radius r2. We find in this case that the solution is given as the 
infinite set of images in two cylinders, 
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As we show in this paper these images are determined completely in terms of q-logarithmic and q-exponential 
functions. Mathematical study of these functions is connected with some applications in the number theory, for 
calculation of the Euler's constant 7 by generalization of a classical formula due to Ramanujan and Vacca |13j . 
irrationality test [12], [13], the Stieltjes transform of a positive discrete measure [11] and construction of the Pade 
approximations [llj . Another class of applications related with physics is quantum groups and their representations 
[16j . The physical systems with such symmetries started from the quantum integrable systems and then extended 
to the several q-deformed physical systems like the quantum linear harmonic oscillator, generalized coherent states 
in quantum optics, composite particle with the Chern-Simons flux - the anyons. It was observed in general that 
physical systems with a fundamental length scale have a symmetry of a quantum group [S] . In nuclear physics the 
deformation parameter is related to the time scale of strong interactions [B], while in solid state physics with the 
lattice spacing [7]. Despite of this particular progress, the direct interpretation of the deformation parameter in 
these cases is sometimes incomplete or even nonexistent. In the present paper we study the classical problem of N 
vortices placed in the annular region between two coaxial cylinders with radii ri < r2 , and we see that the natural 
dimensionless parameter g = r| /rf plays the role of the fundamental length and the infinite set of images produces 
the one dimensional q- lattice. 



2 N Vortices in Annular Domain 

We consider the problem of N point vortices in annular domain D : {ri < \z\ < r2}, where zi, zat are positions 
of vortices with strengths ki, ...,kn respectively. The region is bounded by two concentric circles: Ci : zz = r\ 
and C2 '■ zz = r| 

The complex velocity is given by the Laurent series 

k=l n=0 n=0 

which has to satisfy the boundary conditions at both cylinders, 

[Viz)z + V{z)]\c,^0,k = l,2. (6) 

The conditions ^ imply that no fluid can penetrate any of the circular walls of the domain. If n denotes the 
normal to the boundary, the boundary condition is that normal velocity must be zero, u • v = ii„ = 0. 
To find the unknown coefficients we have to determine the boundary conditions, 

T . Viz)z + Vi2)E = E ^ + E + E ^ + C-C- (7) 

fe=l n=0 n=0 



where C.C. stands for the complex conjugate. Since on boundary Ci: zz = ri^ and \zj^\ > |z|, we can rewrite 
equation ([7]) as foUows; 
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+ C.C. = 0. 



This imphes the following algebraic system 
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0, (n = 0,1,2...) and 5i = 0. 



Since on boundary C2: zz = r2^ and |zfe| < |z|, we can rewrite equation ([7]) as follows; 
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This implies another algebraic system 
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E ^"+2^(^-0^ (n = 0,l,2,...). 

fc=l '"2 '^2 

We have two algebraic systems (jH]) and pT|) . By substracting (fTTj) from ([51), we eliminate a„ 
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li q = r\/r\ is used we find 
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The Taylor series part of ([5]) gives the following, 
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and the Laurent part gives, 



oo , TV oo . 2(n+l) , \ Af oo 

E On+2 _ / 1 -IKk r^' 1 I \^ I 

n=0 k=l \ n=a ■^k ^ I k=l \ n=a 



N / CO . n+1 1 
' 1 IKkZ^ 1 



qn+1 _ I z 



n+1 



N . oo -1 / 9 \ n+1 N . oo 

-iKk M-^ q-l f \ , «Kfc g - 1 



z(q - 1) ^ ff"+i - 1 V zzfe / ^ z(g - 1) ^ - 

fc=l ri=0 ^ ^ ''^ fc=l ^ n=0 ^ 

. oo / 9 \ n+1 N . oo 

Af . . 2 \ ^ ■ 

where [n] = ^fji and Lng{l - x) = - l^^l < 9' > 1- 

3 Logarithmic and q-Exponential Functions 

By analogy with ordinary logarithmic function : < 1, 2; 7^ — 1 
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q-logarithmic function is defined as 
where q-number 
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for any positive integer n. In the limiting case g — > 1 

lim [n] = 71 

lim hnq(\ — x) = ln(l — x) 
Our definition relates with the q-Logarithmic function of Borwein [11] 
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by 

Lng(l - = -Lq{x) 

and with that of Sondow and Zudilin [T5] 
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ln„(l + x) = ^ ^ ^„_^ , \x\<q, q>l 



by 

L7iq(l + x) = (q - 1) ln5(l + x) 

The q-derivative is defined as 

n ^ /(g.T) - /(x) 
= {q-l)x ■ 



Taking q-derivative of a;" and the q-logarithmic function, we get 

Ll^x" = [n]a;""i (28) 

and 



DgLng{l-x) =-j^ (29) 



For our problem, the foUowing representation of q-Logarithmic function [llj . [13j is crucial in the complex 
domain: let q be real, q > 1, then for < \z\ < q the following identity hold 

Ln,{l + z)^Y: ^ Li - ('^ - 1) E (30) 

n=l ^ ' n=l ^ 

The proof is given by the following chain of transformations 
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Next we will need Jackson's q-exponential functions defined as [TT] 
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where [n]! = [l][2]...[n] and q is real. As g — s- 1, both functions reduce to the ordinary exponential. Application of 
q-derivative gives 

DM^) = E<i{x), D,E;{x) = E;iqx) (34) 

The function Eq{z) is entire in z if \q\ > 1, while it has radius of convergence |1 — q\~^ if \q\ < 1. The function 
E*{z) is entire for \q\ < 1 and converges for \z\ < 1 if \q\ > 1. These two functions are related by 

E,{z) = El^^iz) (35) 

and 

E,i~z)E;iz) = 1 (36) 
Particularly important for us is the infinite product representation [llj 
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and the related identity obtained by setting q ^ 1/q and shifting the argument 

n('-^)=T^^'(T^) 
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which is entire for \q\ > 1. 



3.1 Vortex Images and Poles of g-Logarithm 

Substituting equations ITBl) and ([T7]) in ^ we get the following, 
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It can be written in a more compact form in terms of q-derivatives with different basis q: 
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where the real q-number [a]q = (g" — l)/(g— 1), a = a{k) = logq(|zfcp/r^). Parameter a is restricted by < a < f . 
Particular values of this parameter correspond to different positions of a vortex: 1) for a = 1 the vortex is on the 
outer cylinder \zo\ = r2, 2) for a = the vortex is on the inner cylinder \zq\ — ri 3) for a = 1/2 the vortex is at 
the geometric mean distance \zo\ = ^r\r2 (see Section 6 equation ([M)) ) 4) for a = m/n, where m < n are positive 
integers, the vortex is at the generalized mean distance |zo|" — r"^"^r^^. 
Expanding g-log according to ([50)1 we have 
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or 
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Equation for complex velocity has countable infinite number of pole singularities. These singularities can 
be interpreted as vortex images in two cylindrical surfaces. For simplicity let us consider only one vortex at position 



zo, ri < \zq\ < r2- Then the set of images in the cylinder Ci we denote zj 
z^j'j , z j j"* , . . . , where 
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and in the cylinder C2 as 
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Combining together and taking into account alternating signs (the negative for the first image and the positive for 
the next one - the image of the image ) we have two sets of consecutive images 
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This shows that the set of vortex images is completely determined by simple pole singularities of the g-logarithmic 
function. In the above representation by identity r\lq^ = r^/g"^^ we can combine sums so that, we have 
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4 Complex Potential and q-Exponential function 

In this section we derive the complex potential of the flow according to the relation 

V{z)=F\z) i^hZ) 
in terms of Jackson q-exponential function (j32p . To construct F(z) we use new function defined in [llj as 



E, 



where |(?| > 1 and observe 
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Using (|30p we have the relation with q-logarithmic function 

^"^(1-^)- ^.ln/,(z). 
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This expression can be simplified if we use q-derivative of the exponential function ((34| by rescaling the argument 
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By similar arguments we find also 
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If now we apply these formulas to ([39]) we get 
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Finally this implies complex potential in the form 
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The first term in the bracket corresponds to the vortex at position Zk while the second term describes its images. 
All these images completely are determined by zeroes of q-exponential functions. For the zeros of the q-analogue 
of exponential function and asymptotic formulas for varying parameter q see |15) . 



5 Conformal Mapping and Elliptic Functions 

To compare our solution ((6T|) with Johnson & Mc Donald (2004) we rewrite it in terms of elliptic functions. For 

comparison purposes we fix the radius r2 = 1 so that q = ^ = = ^ , where we introduced new parameter 
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q < 1. Then according to ((35)) Eq{z) = Epiz) so that we find complex potential in terms of the second Jackson 
q-exponent 
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Using representation of q-cxponentials as an infinite products (j37p we have 
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The first Jacobi thcta function is defined as an infinite product [17] 
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This theta function is composed from two q-exponcntials 
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Then complex potential becomes 
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By using (|^ we have 
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In terms of coordinates t = —Inz, Tk — —hizk and Tk = — Inzfc, conformally mapping the annulus in the z plane 
to a rectangle in the t plane, finally we find 
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where Fq is a real constant and branch for logarithm is chosen such that ln(— 1) = in, 
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For the stream function, we have 
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This coincides with the result of Johnson & Mc Donald (2004). 



6 Motion of a Point Vortex in Annular Domain 

We use the above formulas to determine the motion of a single vortex in annular domain. Complex velocity at the 
vortex position is determined by, 
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contribution of the vortex itself is excluded. If we take into account that q-harmonic series [14j 

oo 



(76) 



(77) 



converges for g > 1 then at z = zq the first two terms cancel each other and we get the following equation of motion 
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The last equation gives, 
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This implies that the distance of the vortex from the origin is a constant of motion. Then only the argument of 
Zq = |2o|e*'''^*\ is a time dependent function, 

ip{t) =Lut + ipo (80) 



where constant frequency lo dependens on modulus \zq\ 
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So we find that the vortex uniformly rotates around the origin. 
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with frequency depending on the vortex strength, the initial position and geometry of the annular domain. This 
reflects the fact that the motion of vortex results from interaction with an infinite set of its images in the cylinders. 
The frequency (|8ip vanishes when 
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or |zo|"' = rf It means that at the geometric mean distance 
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vortex is at the rest. This equation has simple geometrical meaning that any two intersection points of cylinders 
with a ray (say rie'" and r2e'") are images of each other in a cylinder with radius Y^riri". At this distance angular 
velocity changes the sign and when the vortex approaches the cylinders, lo grows in modulus; w — > K/(2rfe) when 
\zo\ — > rie" « ri(l + e), < e << 1, and lu -K/(2r|e) when \zo\ r2e^'' « r2(l - e), < e << 1. 

Here we like to indicate an intriguing relation of our solution with number theory. The frequency (|8ip is 
combination of two opposite sign frequencies, co = uji + uj2, made from the q-logarithm functions 
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The last representation shows that every frequency is infinite superposition of frequencies coming from every vortex 
image. Moreover for |zo| = t\ the frequency lo\ = H{q), while for |zo| = r2 the frequency lu2 = —H{q) (for simplicity 
we took coefficients equal one), where H(q) is the q-harmonic series. Contribution of N images in the frequency at 
these limiting cases are given by q-harmonic numbers [21j 
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The frequencies (|85|) and ((86)) are compensating each other at the geometric mean distance (|84|. In the annular 
region ri < |zo| < a/^Y^, > \oJ2\ and resulting w > 0, while in the region Y^rir2 < |zo| < ^2, 'j-'i < \^2\ and 
resulting lu < 0. If we consider geometry with parameter q > 2 and the unit strength vortex at the distance 
such that all arguments are non-zero rational, then problem of rationality of the frequency of motion is related 
with problem of rationality of q-logarithms. Starting from early result of Erdos it was proved that the last one is 
irrational [T2]. We can expect influence of this irrationality on the character of multiple vortex dynamics. 



7 The One Vortex Problem in g ^ (X) Limiting Cases 

Usually in applications of q-calculus the limit q ^ 1 corresponds to reduction of the q-elementary functions to the 
standard elementary functions. However in our problem this limit corresponds to ri — r2 and the region reduces 
to the circle. More interesting is the limit when q oo. To study this limit we need corresponding limits of 
q-elementary functions. When the q-logarithm is expanded we get the following 
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so that 







z2 

I TT7T = 1 + Z A 

0"("-l)/2 Q 
n=0 ^ ^ 



lim E„{z) = 1 + z. 

q — >oo 



Applying the q-dcrivative 



we have identities 



E,iqz)-Eg{z) 



E. 



{q-l)z 

qz 

Ea 



DgEgiz) ^ Eg{z), 



In the hmit q ^ 00 they imply that 



lim Eg 



= il+z)Eg 

Eqjz) 
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= 1 



lim E„ 

q — ^00 



qz 

q-l 



1 + z. 



(93) 
(94) 

(95) 

(96) 

(97) 
(98) 

(99) 
(100) 



Since q — r^/rl, for the limit g — > 00 we have two different geometrical cases: 

1) when ri = constant and r| = qr\ 00, the outer cylinder grows unlimited so that we have just the one 
cylinder problem. Complex velocity ([55]) by replacement — qr\, taking the limit q ^ 00 and using (|89p . (PT|) 
gives finite result 

IK in r\/zQ 



V{z)^ 



(101) 



z - zo z z-rl/zo 

which is exactly the result of the circle theorem ^ . For angular velocity we have using (|8T|) and the limits ([M)) , (191]) , 

^2 

For the complex potential we have 

F{z) — inln 

Applying the limits ((99)l . (|100p we obtain two images as required by the circle theorem 

F{z) = in 



{z - 20)- 



E„ 



E„ 



'1 \(l-q)zo y ^9 \(l-q)z 



qrf 



(-3) 


+ In z 







(103) 



(104) 



2) when r2 = constant and = r|/g — > 0, the inner cylinder decreases until the thin string and then 
disappears. In this limiting case in complex velocity ([39]) we replace rf = r^/g, and use formulas ([89 ]) . ((9T|) so that 
we get expression ([3]) . For angular velocity of one vortex from (|8T1) by limits (j89|) , ([91]) we have 

Applying above limits ((99| . (|100p to the complex potential 



F{z) — IK In 



{z - zo)- 



E„ 



'1 \{l-q)zo ; ^1 \ 



En 



E„ 



1 V (l-q)ri ) ^1 



\ (i-g)2zo ) 



we obtain just one image as expected from ([3]) 



F{z) = in 



ln(z - 20) ~ In U - ^ + In 



zq 



zo 



(106) 



(107) 



8 Conclusions 



We have shown that the infinite set of images for a vortex in annular domain between two coaxial cylinders is 
given completely in terms of q-logarithmic function for complex velocity and in terms of the Jackson q-exponential 
function for the complex potential. Recent results on Pade approximation for the q-elementary functions [11] could 
be efficient approximation in the vortex image problem, restricting number of images. For example in paper [13] 
in order to compute the q-logarithm, with reference on P. Sebah, the next formula proposed 



N 

Zn,(l + z) = z^^^+r^(z) (108) 

n—1 



where 



n=l ^ 

with N being any positive integer. Then application of this formula to our problem gives N-vortex images approx- 
imation which could be useful in applications. 

By conformal mapping of the annular region to the rectangular one and composing q-exponents, the stream 
function has been represented in terms of the Jacobi elliptic function and it coincides exactly with result of 
Johnson, Mc Donald. Since annular region can be conformally mapped to exterior of two cylinders in the plane, 
our solution provides also solution of the last problem in terms of q-elementary functions. Results of calculations 
are in preparation [10 . 

Finally we like to note that elementary relation between the hydrodynamical problem and q-calculus considered 
in this paper could be applied to several physical situations with the same geometry, as the electrostatic problem 
or the problem of anyons [8 . The image picture allow us to construct Green function in the domain in terms of 
q-elementary functions and apply it to other problems like the vortex dynamics. 

9 Acknowledgements 

This work was partially supported by TUBITAK under the Grant No. 106T447 and by Izmir Institute of Technol- 
ogy- 
References 

[1] Milne-Thomson L. M., Theoretical Hydrodynamics, London, Macmillan, 1968 
[2] Power G., A dielectric cylinder. Math. Gaz. 37 (1953) 220 

[3] Palaniappan D., Classical image treatment od a geometry composed of a circular conductor partially merged 
in a dielectric cylinder and related problems in electrostatics, J. Phys. A: Math. Gen. 38 (2005) 6253-6269 

[4] Yilmaz, O., An Iterative Procedure for the Diffraction of Water Waves by Multiple Cylinders, Ocean Engi- 
neering, 31, 1437-46, (2004). 

[5] Bonechi F., Celeghini E., Giachetti R., Sorace E. and Tarlini M., Quantum Galilei Group as symmetry of 
magnons Phys. Rev. B46 (1992) 5727 

[6] Celeghini E., Giachetti R., Sorace E. and Tarlini M., Quantum groups of motion and rotational spectra of 
heavy nuclei, Phys. Lett. B (1992) 

[7] Bonechi F., Celeghini E., Giachetti R., Sorace E. and Tarlini M., Inhomogencous Quantum Groups as sym- 
metry of phonons, Phys. Rev. Lett. 68 (1992) 3718 

[8] Pashaev O. K. , Gurkan Z. N., Abelian Chern- Simons Vortices and Holomorphic Burger's Hierarchy, Theor. 
Math. Phys.l52{l) (2007) 1017-1029 . 

[9] Johnson E. R., McDonald N. R., The Motion of a Vortex Near Two Circular Cylinders, Proc. R. Soc. London 
A, 460, 939-54, (2004). 



[10] Cagatay F., Gurkan Z. N., Pashaev O.K., Yilmaz O., Vortex Images in Two Cylinders as Poles of q-Logarithmic 
Function, to be published. 



[11 
[12 

[13 

[14 
[15 

[16; 

[17: 

[18 
[19 

[20; 

[21 



Borwein P. B., Pade Approximation for the q-Elementary Functions, Constr. Approx. 4 (1988) 391-402 
Borwein P. B., Journal of Number Theory, 37 (1991) 253-259 

Sondow J., Zudilin W., Euler's Constant, q-Logarithms, and Formulas of Ramanujan and Gosper, The Ra- 
manujan Journal, 12 N2 (2006) 225-244 

Zudilin W., Approximations to q-logarithms and q-dilogarithms, with applications to q-zeta values, (2004) 

C A Nelson and M G Gartley, On the zeros of the q-analogue exponential function, J. Phys. A:Math Gen 27 
(1994) 3857-3881 

A. Klymik and K. Schmudgcn, Quantum Groups and their representations. Springer, 1997 

E. T. Whittaker, G. N. Watson, A Course of Modern Analysis, fourth edition, Cambridge Univ Press, London, 
1963 

R. H. Levy, Diocotron instability in a cylindrical geometry. Physics of Fluids, 8 (1965) 1288 

M. Amoretti, D. Durkin, J. Fajans, R. Pozzoli and M. Rome, Non-Neutral Plasma Physics IV, edited by F. 
Anderedd et al., AIP, 2002 

Gibbons G.W., Vortex motion in the neighbourhood of a cosmic string, Proc. R. Soc. Lond. A 425 (1989) 
431-439 



Shi Ling-Ling, Pan Hao, A q-analogue of Wolstenholme's harmonic series congruence, |math.NT/0507495| 



